Tailoring optical nonlinearities via the Purcell effect 
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We predict that the effective nonlinear optical susceptibility can be tailored using the Purcell 
effect. While this is a general physical principle that applies to a wide variety of nonlinearities, 
we specifically investigate the Kerr nonlinearity. We show theoretically that using the Purcell 
effect for frequencies close to an atomic resonance can substantially influence the resultant Kerr 
nonlinearity for light of all (even highly detuned) frequencies. For example, in realistic physical 
systems, enhancement of the Kerr coefficient by one to two orders of magnitude could be achieved. 
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Optical nonlinearities have fascinated physicists for 
many decades because of the variety of intriguing phe- 
nomena that they display, such as frequency mixing, 
supercontinuum generation, and optical solitons [TJ 
Moreover, they enable numerous important applications 
such as higher-harmonic generation and optical signal 
processing [3J[31|3]. On a different note, the Purcell ef- 
fect has given rise to an entire field based on studying 
how complex dielectric environments can strongly en- 
hance or suppress spontaneous emission from a dipole 
source [5J El |5J W\ ■ In this Letter, we demonstrate that 
the Purcell effect can also be used to tailor the effective 
nonlinear optical susceptibility. While this is a general 
physical principle that applies to a wide variety of nonlin- 
earities, we specifically investigate the Kerr nonlinearity, 
in which the refractive index is shifted by an amount 
proportional to intensity. This effect occurs in most ma- 
terials, modeled here as originating from the presence 
of a collection of two-level systems. We show theoret- 
ically that using the Purcell effect for frequencies close 
to an atomic resonance can substantially influence the 
resultant Kerr nonlinearity for light of all (even highly 
detuned) frequencies. 

In hindsight , the modification of nonlinearities through 
the Purcell effect could be expected intuitively: opti- 
cal nonlinearities are caused by atomic resonances, hence 
varying their strengths should influence the strengths of 
nonlinearities as well. Nevertheless, to the best of our 
knowledge, this interesting phenomenon has not thus far 
been described in the literature. Moreover, as we show 
below, it displays some unexpected properties. For ex- 
ample, while increasing spontaneous emission strength- 
ens the resonance by enhancing the interaction with the 
optical field, it actually makes the optical nonlinearity 
weaker. Furthermore, phase damping (e.g., through elas- 
tic scattering of phonons), which is detrimental to most 
optical processes, plays an essential role in this scheme, 
because in its absence, these effects disappear. 

A simple, generic model displaying Kerr nonlinearity 
is a two-level system. Its susceptibility has been cal- 
culated to all orders in both perturbative and steady 
state limits [2]. However, this derivation is based on a 



phenomenological model of decay observed in a homoge- 
neous medium, and does not necessarily apply to systems 
in which the density of states is strongly modified, such 
as a cavity or a photonic crystal bandgap. Following an 
approach similar to Ref. [TOl the validity of this expres- 
sion can be established from a more fundamental point of 
view. Start by considering a collection of N two-level sys- 
tems per unit volume in a photonic crystal cavity, whose 
levels are labeled a and b. The corresponding Hamilto- 
nian is given by the sum of the self-energy and interac- 
tion terms (Hq and V(t), respectively). Using the electric 
dipole approximation, one obtains: 

H = H +V(t) = h [uj a a aa + u b a bb + fl(t)a ab + Q* (t)a ba ] , 

(1) 



where <7,- 



c\cj is the operator that transforms the 
fermionic state j to the fermionic state i, £l(t) = —fx ■ 
E{t)/h is the Rabi amplitude of the applied field as a 
function of time, and the scalar dipole moment p, is de- 
fined in terms of its projection along the applied field 
E(t) . In general, if this system is weakly coupled to the 
environmental degrees of freedom, then the timescale for 
the observable dynamics of the system is less than the 
timescale of the "memory" of the environment. In this 
case, information sent into the environment is irretriev- 
ably lost - this is known as the Markovian approxima- 
tion [TT] . The dynamics of this system can then be mod- 
eled by the Lindbladian C, which is a superoperator de- 
fined by p = C [p] . In general, one obtains the following 
master equation from the Lindbladian: 



= -(i/ft) [H,p] + J2 
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(2) 

Using the only two quantum jump operators that are al- 
lowed in this system on physical grounds - L\ = o ab j ^/T\ 
and L 2 = o- bb ^/"f pheiBB [llj - one can obtain the following 
dynamical equations: 



dp ba 
dt 



-(ioj ba + T 2 1 )p ba + in(t)(p bb - p aa ) (3) 
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bb 



dt 



— 2i [n(t)p ab ~n*(t) Pba ]. 

(4) 

where uj ba = uj b — ui a , T\ is the rate of population loss 
of the upper level, and T 2 — (l/2)T^ 1 + 7 phaao is the 
rate of polarization loss for the off-diagonal matrix el- 
ements. The prediction of exponential decay via spon- 
taneous emission is known as the Wigner-Weisskopf ap- 
proximation [12]. Although it has been shown that the 
atomic population can display unusual oscillatory be- 
havior in the immediate vicinity of the photonic band 
edge [T3J [U], theoretical [TU] and experimental consid- 
erations |15l 116] show that this approximation is fine for 
resonant frequencies well inside the photonic bandgap. In 
the rest of this manuscript, this is assumed to be the case. 
Next, one can make the rotating wave approximation for 
Eqs. (J3j) and Q, and then solve for the steady state. If 
the polarization is defined by P = N/j,(pb a + p a b) = xE, 
where x is the total susceptibility to all orders, one ob- 
tains the following well-known expression for the suscep- 
tibility 0E!]: 



X = 



Np 2 (io~io ba -t/T 2 )Ti/H 
l + (w- coba^Ti + (4/tP)n*\E\*T l T 2 ■ 



(5) 



In general, equation ([5| may be expanded in powers of 
the electric field squared. Of particular interest is the 
Kerr susceptibility, also in Ref. |2] 
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where A = to — oj ba is the detuning of the incoming wave 
from the electronic resonance frequency. For large detun- 
ings AT2 3> 1, one obtains the approximation that: 



Re X (3) « -Nn A 
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Of course, there are many types of materials to which 
a simple model of nonintcracting two-level systems does 
not apply. However, it has been shown that some semi- 
conductors such as InSb (a III-V direct bandgap mate- 
rial) can be treated as a collection of independent two- 
level systems with energies given by the conduction and 
valence bands, and yield reasonable agreement with ex- 
periment |17j . If the parameter A is defined in terms of 
the bandgap energy such that A = ujq — u>, then one can 
look at the regime AT 2 ^> 1 studied above, and obtain 
the following equation: 
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where P is a matrix element discussed in Ref. |T7j loq 
is the direct bandgap energy of the system, and m r is 
the reduced effective mass of the exciton. This equation 
displays the same scaling with lifetimes as Eq. Q, so 



the considerations that follow should also apply for such 
semiconductors . 

Now, consider the effects of changing the spontaneous 
emission properties for systems modeled by Eqs. ([T]) or 
When spontaneous emission is suppressed, as in 
the photonic bandgap of a photonic crystal, Ti will be- 
come large while T 2 remains finite, thus enhancing x^ 
by up to one or more orders of magnitude (for materials 
with the correct properties). For large detunings (where 
AT 2!Vac > 1), we expect that x (3) will scale as Ti/T 2 . 
The enhancement of the real part of x^ 3 ' is defined to 
be r, = Re^leu/Rexll where X £L» is the nonlin- 

ear susceptibility in the presence of the Purcell effect, 

(3) 

while xLm is the nonlinear susceptibility in a homoge- 
neous medium. Since Tf 1 = T rad + T nr , the maximum 
enhancement is predicted to be roughly: 



^r,,,. + 7 Phas 



r rad + r nr 



11 ~ r 1)hom r 2)Purcoll ~ i(r nr + r rad ) + 7phaac 

(9) 

where r rad is the radiative decay rate in vacuum. Since 
the Purcell effect increases the amplitude of x^\ one 
might also expect it to increase the amplitude of x^', 
however, according to Eq. the opposite is true. This 
can be understood by noting that Purcell enhancement 
decreases the allowed virtual lifetime, and thus, the like- 
lihood of nonlinear processes to occur |18j . Moreover, 
since the Purcell factor [S] is calculated by only consid- 
ering the photon modes [?,, one would not necessarily 
expect phase damping effects to also play a role. How- 
ever, the results of Eq. (|9| show the contrary to be true, 
and can be explained as follows: when phase damping 
is large, the polarization will decay quickly, thus giving 
rise to a small average polarization. However, as phase 
damping is lessened, polarization decay slows down and 
allows the average polarization to rise. In the limit that 
phase damping is controlled exclusively by the sponta- 
neous emission rate, the two competing effects will can- 
cel, and the nonlinearity will revert to its normal value. 

Furthermore, the presence of large phase damping ef- 
fects makes T 2 effectively constant, which means that 
suppression of spontaneous emission (caused by the ab- 
sence of photonic states at appropriate energies [5]) can 
enhance Kerr nonlinearities by one or more orders of 
magnitude, while enhancement of spontaneous emission 
via the Purcell effect [7J |5] can suppress these nonlin- 
earities. For the case where Purcell enhancement takes 
place, Ti decreases while T 2 may not change as rapidly, 
due to the constant contribution of phase damping ef- 
fects. This applies in the regime where T\ 3> 7^^ e . 
Otherwise, for sufficiently small Ti, T 2 will scale in the 
same way and x*' 3 '' will remain approximately constant 
for large detunings, where AT 2 3> 1. This opens up 
the possibility of suppressing nonlinearities in photonic 
crystals (to a certain degree). For processes such as 
four-wave mixing or cross phase modulation, x^ will 
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generally involve a detuning term and will differ from 
Eq. dfjl. It is also interesting to note that this enhance- 
ment scheme will generally not increase non-linear losses, 
which are a very important consideration in all-optical 
signal processing. If the nonlinear switching figure of 
merit £ is defined by £ = Re%( 3 V(AImx^) P^l - tnen 
Cpurocii/^vaouum = T 2 , purcc „ / T 2 .vacuum > 1, for all cases of 
suppressed spontaneous emission. 

The general principal described thus far should apply 
for any medium where the local density of states is sub- 
stantially modified. In what follows, we show how this 
effect would manifest itself in one such exemplary system: 
a photonic crystal. This example serves as an illustration 
as to how strong nonlinearity suppression / enhancement 
effects could be achieved in realistic physical systems. It 
consists of a 2D triangular lattice of air holes in dielectric 
(e = 13), with a two-level system placed in the middle, 
as in Fig. [T] 

Note that the vast majority of photonic crystal liter- 
ature is generally focused on modification of dispersion 
relations at the frequency of the light that is sent in as a 
probe. By contrast, in the current work, it is only essen- 
tial to modify the dispersion relation for the frequencies 
close to the atomic resonances; the dispersion at the fre- 
quency of the light sent in as a probe can remain quite 
ordinary. 




FIG. 1: (Color) A 2D triangular lattice of air holes in dielec- 
tric (e = 13). On top of the dielectric structure in grey, the 
E z field is plotted, with positive values in red, and negative 
values in blue. A small region of nonlinear material is placed 
exactly in the center of the structure. This material may be, 
for example, either two-level atoms, quantum wells, or some 
semiconductors such as InSb. 

First, consider the magnitude of the enhancement 
or suppression of spontaneous emission in this system. 
Clearly, since there are several periods of high contrast 
dielectric, two effects are to be expected. First, there 
will be a substantial but incomplete suppression of emis- 
sion inside the bandgap. Second, there will be an en- 
hancement of spontaneous emission outside the bandgap 
(since the density of states is shifted to the frequencies 
surrounding the bandgap). For an atom polarized in the 
direction out of the 2-D plane, only the TM polarization 
need be considered. 



We numerically obtain the enhancement of sponta- 
neous emission by performing two time-domain simula- 
tions in Meep [20], a finite difference time-domain code 
which solves Maxwell's equations exactly with no approx- 
imations, apart from discretization (which can be sys- 
tematically reduced) [2T]. First, we calculate the spon- 
taneous emission of a dipole placed in the middle of the 
photonic crystal structure illustrated in Fig. [TJ then di- 
vide by the spontaneous emission rate observed in vac- 
uum. The resulting values of T\ and T 2 are calculated 
numerically, and Eq. ([6]), in conjunction with the defini- 
tion of the enhancement factor 77, is used to plot Fig. [3] 
The results are plotted in Fig. [2] 
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FIG. 2: (Color) Relative enhancement of the TM local den- 
sity of states for Fig. [TJ as measured in the time-domain sim- 
ulation rate of emission, F, normalized by the emission rate 
in vacuum, r o 

A GaAs-AlGaAs single quantum well can lie in the 
interesting regime discussed above, where the radiative 
loss rate r rad dominates the non-radiative loss rate r„ 
as well as the overall loss rate of the quantum well, for 
certain temperatures |22j . Equation ^ implies that one 
can see substantial enhancement of the Kerr coefficient 
in that regime. 

At a temperature of about 200 K, r nr w 0.ir rad [22]. 
Although experimental measurements for 7 phaBC are un- 
available to the authors, the presence of a substantial 
phonon bath at that temperature leads one to expect a 
fairly large value, which may be conservatively estimated 
by 10r rad . These results are displayed in Fig. [3ja). Note 
that enhancement is primarily observed inside the pho- 
tonic bandgap (cf. Fig. [2]). We observe an enhancement 
in the real part of the Kerr coefficient up to a factor of 
12, close to the predicted maximum enhancement factor 
of 10.48 in the regime of large detunings (AT2 3> 1). 

Also, at a temperature of about 225 K, r„ r ss r rad [22J, 
and again we take 7 pha3C = 10r rad . These results are 
displayed in Fig. |3jb). In this case, we observe an en- 
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FIG. 3: (Color) Contour plot of Kerr enhancement n = 

R e Xt>urceii/ R- e XhoL as a function of probe (uj ph ) and elec- 
tronic transition (a> e iec) frequencies, for a single quantum well 
of GaAs-AlGaAs, (a) at T = 200 K, with 0.l7 phasc = 10r nr = 
r rad , and (b) at T = 225 K, with 0.l7 pha5C = P„ = r rad . 



hancement up to a factor of 2.5, close to the predicted 
maximum enhancement factor of 1.91 in the regime of 
large detunings (AT 2 S> 1). 

Finally, we note that close to room temperature (285 
K), the system in Rcf. 22 displays F nc sa 10r rad , which 
is predicted to yield a maximum enhancement factor of 
1.06. Since this number is fairly negligible, it illustrates 
that this approach has little impact when non-radiative 
losses dominate the decay of the electronic system. 

On the other hand, some recent work has demon- 
strated that a single quantum dot can demonstrate pre- 
dominantly radiative decay in vacuum even at room 
temperature, e.g., single CdSe/ZnS core-shell nanocrys- 
tals with a peak emission wavelength of 560 nm, with 
T ra d ~ 39r„ r . Even bulk samples of similar nanocrys- 
tals have been shown to yield a significant radiative decay 
component, corresponding to r rad « T nr |24j . Thus, we 



predict that with strong suppression of radiative decay, 
nonlinear enhancement of a factor of two, or more, could 
be observed at room temperature. 

We now discuss the implications of this effect on pre- 
vious work describing nonlinearities in photonic crystals, 
such as Ref. [25j and the references therein. Most past ex- 
periments should not have observed this effect, because 
they were designed with photonic bandgaps at optical 
frequencies significantly smaller than the frequencies of 
the electronic resonances generating the nonlinearities, 
in order to operate in a low-loss regime. Furthermore, in 
most materials, non-radiative decays will dominate ra- 
diative decays at room temperature. Finally, all the pre- 
vious analyses are still valid as long as one considers the 
input parameters to be effective nonlinear susceptibili- 
ties, which come from natural nonlinear susceptibilities 
modified in the way described by this paper. 

In conclusion, we have shown that the Purcell effect 
can be used to tailor optical nonlinearities. This principle 
manifests itself in an exemplary two-level system embed- 
ded in a photonic crystal; for realistic physical param- 
eters, enhancement of Kerr nonlinearities by more than 
an order of magnitude is predicted. The described phe- 
nomena is caused by modifications of the local density of 
states near the resonant frequency. Thus, this treatment 
can easily be applied to analyze the Kerr nonlinearities 
of two-level systems in almost any geometrical structure 
in which the Purcell effect is substantial (e.g., photonic 
crystal fibers [55] , optical cavities). It also presents a 
reliable model for a variety of materials, such as quan- 
tum dots, atoms, and certain semiconductors. Future in- 
vestigations will involve extending the formalism in this 
manuscript to other material systems. 
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